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Abstract. We outline the scientific motivation behind a search for gravitational waves 
associated with short gamma ray bursts detected by the Interplanetary Network (IPN) during 
LIGO's fifth science run and Virgo's first science run. The IPN localisation of short gamma 
ray bursts is limited to extended error boxes of different shapes and sizes and a search on these 
error boxes poses a series of challenges for data analysis. We will discuss these challenges and 
outline the methods to optimise the search over these error boxes. 



1. Introduction 

Short hard gamma ray bursts (short GRBs) are believed to be produced by mergers of either 
two neutron stars or a neutron star and a stellar mass black hole [HEIE]. These events are ideal 
sources for strong gravitational wave emission [H [5] . If an observation of both gamma rays and 
gravitational waves (GW) originating from the same event could be achieved, it will increase 
confidence and allow for better science output. It is thus of great importance to constantly 
monitor and record GRBs to allow a GW search to be performed around the burst times. 
Systematic analyses of GW data around short GRB times have been done in the past and the 
most recent publications from the LIGO- Virgo group contain results from the Swift-observed 
GRBs during LIGO's fifth science run (S5) and Virgo's first science run (VSRl) [HI [7]. This 
paper reports the methodology and motivation for a proposed GW search around the times of 20 
additional GRB during S5/VSR1. These bursts were observed by the Interplanetary Network 
(IPN) [HI |9], a group of satellites orbiting the Earth and Mars and operating, among other 
equipment on board, gamma ray detectors. These bursts were detected between 2006 and 2007 
and have been localized, in both time and sky location, such that a targeted GW search is 
possible. 

The Interplanetary Network [SJ [TU] employs several space missions and synthesizes data 
obtained from the detection of the same burst by different spacecraft equipped with gamma 
ray detectors. The IPN has been operating for three successive generations; presently the 
third IPN (IPN3) began its operation in November 1990. Currently the spacecraft gathering 
data are Konus-WIND, Suzaku, INTEGRAL, RHESSI, Swift, Fermi/GBM (in Earth orbit), 
MESSENGER (in Mercury orbit) and Mars Odyssey (in Mars orbit) [9|. When the duty cycles 
and effective fields of view of all the missions in the network are considered, the IPN is an 
all-time, isotropic GRB monitor. 



In this paper, we first review the motivation for performing a joint GW-GRB search as an 
overview and then provide details of the IPN short GRBs that were identified during the LIGO- 
Virgo S5-VSR1 science runs. We discuss the necessary changes to the current ongoing analysis 
(for the S6/VSR2-3 science runs Swift and Fermi-observed GRBs, [Ulllll) that are needed to 
implement a search on GRBs identified by the IPN network which may be less well localized in 
both sky position and time than the corresponding bursts identified by the Swift satellite and 
used in previous analyses [Gj [7] . 

2. Search for GW associated w^ith short GRB 

Gamma ray bursts fall into two commonly accepted groups depending on their duration 
and spectral hardness [U |13]. Characteristic duration is quantitatively expressed by the Tgo 
parameter, defined as the time interval over which 90% of the total background-subtracted 
counts are observed (the total counts emitted by the source, found from the counts in the 
source region minus the contribution from the background), with the interval starting when 5% 
of the total counts have been observed [14^. The majority of bursts, with softer spectrum and 
duration Tgo > 2s are called long GRBs. They are thought to be produced during a core-collapse 
supernova event. Those with a harder spectrum and Tgo < 2s are called short GRBs. 

Short GRBs (Tgo < 2s) are widely believed to be mergers of either neutron star-neutron 
star or neutron star-black hole binaries [T]. Such compact binary coalescences are believed to 
generate strong gravitational waves in the sensitive frequency band of Earth-based gravitational- 
wave detectors |15l 116) . Although binary coalescence is the favoured progenitor model for short 
GRBs, this has not yet been confirmed by means of a direct observation. Consequently, the 
detection of gravitational waves associated with a short GRB would provide evidence that the 
progenitor is indeed a coalescing compact binary and also provide information on the parameters 
of the binary. 

The coalescence of a compact stellar-mass binary (either two neutron stars or a neutron star 
and a black hole companion) is the endpoint of its ~Gyr life revolving around the common centre 
of mass while constantly losing energy and angular momentum through emission of gravitational 
radiation. The binary merger will lead to the formation of either a transient hyper massive 
highly magnetized neutron star with a lifetime of a few ms \n\ [18] that will collapse to a rapidly 
spinning black hole or straight to the formation of a black hole [21 [19] . In the favoured model 
of short GRBs [ 17[ [SJ [20l I21| . the gamma-ray emission is contingent on the formation of a 
massive (and highly magnetized) torus around the final black hole. The matter in the torus is 
accelerated to relativistic velocities leading to the formation of a coUimated jet of electromagnetic 
radiation along the axis of the former binary total angular momentum. Differences in velocities 
between layers in the jet account for the prompt gamma-ray emission. Full relativistic and 
magnetohydrodynamic numerical simulations have shown that the time difference between the 
binary merger and the jet formation can be between a few milliseconds up to a few seconds 
[21 [131201 [21]. 

Supposing gravitational waves produced just before merger travel at the same speed as the 
gamma rays, the speed of light in vacuum, and suppose the observer is situated within the cone 
of the jet, one would expect to observe the GW within a few seconds prior to the arrival of the 
gamma rays. As the GWs emitted during the pre-merger inspiral dominate the signal-to-noise- 
ratio (SNR) observable in current GW detector data and they can be accurately modelled using 
post-Newtonian approximations [221 [23] , our search will be aimed at this phase only. 

3. The short GRB sample to be analysed 

We gathered the sample of GRBs to be analysed using data provided in the IPN online table 
available publicly at [9j and by cross-checking this with NASA's HEASARC online table at [M] . 
Since these GRBs are not reported in any of NASA's Gamma Ray Burst Circulars (GCN), a 



manual check on each burst had to be performed in order to confirm its characteristics, involving 
checking several IPN satellite homepages e.g. Suzaku at INTEGRAL at [26], Swift at [27] 
and Konus-WIND at pS]. 

The details of the sky position were obtained by manually processing raw data from the 
IPN satellites for every GRB. In order to perform a search of the GW data for a given GRB, 
the sky position and time of the GRB must be determined. This information is obtained from 
knowledge of: the IPN satellites that detected the burst together with their absolute and relative 
sky positions (information needed for constructing the GRB error boxes), locations of all the 
spacecraft (used to obtain the blocking constraints to reduce the size of the GRB error boxes), 
the burst time of arrival at the satellites and at Earth (used to determine the time interval 
on which we will perform the GW search) and its error. Constructing the GRB error boxes, 
determining their sizes and choosing the right GW data for analysis are entirely contingent on 
these pieces of information. 



Figure 1. The IPN schematics: triangulating the position of a GRB using three IPN spacecraft 
(SI, S2 and S3). Using three satellites we obtain two IPN annuli that intersect to form two error 
boxes. Reference P 



3.1. The IPN triangulation mechanism 

The principle on which the IPN triangulation method is based uses the arrival time of the same 
burst at different spacecraft to determine the source sky location. Figure [l] illustrates how an 
IPN triangulation works. SI, S2 and S3 denote three spacecraft detecting the same GRB and 9 
is angle between the burst direction and the baseline between SI and S2. Then, the burst will 
be detected by S2 at a time 5T seconds earlier than SI 



where D12 is the distance between SI and S2, and c the speed of light. Since D12 is known and 
5T is measured, 9 is estimated. The solution to equation (1) is represented by a ring, or an 
annulus, whose width depends on the timing uncertainties {a{6T)) and on the separation Di2. 
The farther apart the detectors, the more precise the localisation. The number of independent 
couples of detectors (and, therefore, the number of independent annuli) is two for the case of 
three spacecraft; thus, the burst direction must be inside one of the two intersection regions 
of the annuli. This intersection region is called the IPN error box of the GRB. Depending on 




ITHE " TRIANGULATION" METHOD] 



cos{9) = C6T/D12 



(1) 



the location of the IPN spacecraft and their timing errors this error box may vary in size from 
fractions of, to hundreds of square degrees. The annulus width is obtained by propagating the 
uncertainty on the time delay 5T. Thus, from equation (1) it follows that 

Equation ^ gives the uncertainty a{9) in the angle 9 expressed in radians. The uncertainty 
on Di2 has been neglected, as the main contribution comes from the timing uncertainties. One 
has to take into account not only the time resolution of each of the detectors, but also the 
difficulty of comparing different light curves, often derived from different energy bands. For 
example, when D12 spans the typical range: few 10^ - few 10^ light seconds, then from equation 
([2]) it comes out that a minimum time resolution of the order of 10~^ — 10~^ s is required, in 
order to have a{9) ~ few arcminutes or less in order to obtain a precise localisation. 




Figure 2. Error box construction: red lines correspond to the Konus-WIND ecliptic latitude 
bands, blue lines are the 3-0" IPNl-2 annulus and magenta lines are the 3-0" IPNl-3 annulus 
obtained by triangulation (in this case we have three IPN spacecraft observing the same burst). 
The error boxes are the solid regions bounded by the intersections of the IPN annuli and ecliptic 
bands. 



3.2. The IPN GRB error box construction 

The error boxes for the IPN GRBs are constructed from the intersection of the triangulated 
3-0" IPN annuli and different field-of-view blocking constraints (if present). See figure [2] for an 
illustration example. When the constraints involve Konus-WIND ecliptic latitude bands, the 
error box that will be kept is located between the ecliptic bands. The reason for this is that 



the Konus-WIND spacecraft consists of two detectors, one facing the north echptic pole, and 
the other facing the south echptic pole. By comparing the count rates on these two detectors, 
the Konus team obtains an estimate of the ecliptic latitude of the burst. Typically, the band 
is 20-30 degrees wide, and it is good to a 90-95% confidence; systematic and statistical errors 
usually prevent it from being much more accurate than this value. In those cases where we have 
a single triangulation annulus, plus an ecliptic latitude band, the result is often two long, narrow 
error boxes, where the annulus intersects the band. Other intersections are possible, too, such 
as a single long, narrow error box if the IPN annulus grazes the ecliptic latitude band. Also, 
anything that overlaps a region where a planet blocks the view will be discarded. 

3.3. Determining the GRB time of arrival 

The Earth crossing time, also referred to as the time of arrival of the burst, is the time when 
the gamma ray signal would cross the centre of the Earth and is the reference time to be 
considered when constructing the time search window for gravitational waves. This time can 
only be estimated based on the burst arrival times at the different IPN satellites and based on 
their positions with respect to Earth. This way of choosing the burst time is prone to two types 
of uncertainties: the first is simply the fact that we may not have the time at Earth but at a 
satellite that is separated from Earth by a certain distance; the error is directly proportional to 
the distance to Earth where the closest IPN satellite is located at the time of the burst. The 
best estimate comes from any satellite that is not interplanetary (i.e. not on orbit around Mars 
or at a distant point from Earth). These "close" satellites are usually within a fraction of a 
light second distance from Earth. These uncertainties are small, less than 1 second for satellites 
orbiting Earth but may be of up to 5 seconds or more for interplanetary satellites. The second 
kind of uncertainty comes from the fact that the IPN consists of nine spacecraft with different 
energy ranges and different spectral sensitivities. So it is possible, in an extreme case, to have 
one spacecraft trigger on a GRB precursor, while the others trigger on the main burst, resulting 
in a time difference. This way, the trigger times can differ by tens of seconds or more. For 
short GRBs this effect is minimised due to the hard nature of their spectra and consequently 
reduces to under one second imprecisions. Altogether, time imprecisions are not more than 5s 
for the burst sample we will be analysing. The few GRBs that have a timing imprecision greater 
than 1 second are localised by distant satellites, usually MESSENGER/Mars Odyssey and/or 
Konus-WIND. 

3.4-. Gravitational wave detectors and data availability 

Our aim is to perform a search for GW associated with the well or fairly well localised short 
GRBs detected by the IPN during LIGO's S5 run that lasted from 4 November 2005 to 30 
September 2007, and Virgo's VSRl that lasted from 18 May 2007 to 30 September 2007. The 
analysis will make use of data from four operational GW interferometers: the 4km and 2km 
co-located LIGO detectors at Hanford, WA (HI and 112), the 4km detector at Livingston, LA 
(LI) and the 4km detector at Cascina, Italy (VI) [29l [301 [31] . The search will be using a method 
that coherently combines data from multiple operational GW detectors described in [12]; the 
method is being used in GW-GRB searches for S6/VSR2-3 |llj and is proven to be performing 
better than the method used for the S5/VSR1 search [121 [3 [7j; we require that all GRBs have 
data from at least two GW detectors. In order to perform the search, we require approximately 
forty minutes of data around the time of the GRB. The search pipeline identifies a foreground 
time representing the time interval around the actual burst when the signal is most likely to have 
been detected by the GW detectors. For the IPN GRBs that have a burst time of arrival error 
less than a second, based on the delay between the time of the arrival of the gamma ray and of 
the GW (described above) an interval of 5 seconds prior to the GRB time and 1 second following 
it will be used as foreground. This time interval will account for any timing imprecisions and 



has been previously used in the S5/VSR1 search for GW associated with the Swift GRBs [H]. 
For the IPN GRBs that have a time of arrival error larger than 1 second, the foreground will be 
extended on either side to account for this error. The data surrounding the time of the GRB 
are used for background estimation, in order to assess the data quality in the detectors around 
the time of the GRB. 

The LIGO- Virgo detector network |29| [30l [5T] is sensitive to a large fraction of the sky, albeit 
with relatively poor localization capability (on the order of tens or hundreds of square degrees) 
|32j . In much the same way as the IPN network, the GW network of detectors can reconstruct 
a sky position primarily through triangulation. A GW signal with an SNR large enough to be 
detected will have its sky location resolved by timing its arrival at the different GW detectors. 
The approximate timing resolution between two detectors in the network is ~0.5ms, giving a 
best angular resolution of around 2°. When performing a search for a GW signal associated 
with a GRB, the data from the various detectors in the network is coherently added with the 
appropriate time delays between detectors corresponding to a given sky location. Thus, if the 
IPN error box spans a large region of the sky, it is necessary to search several different sky 
positions for a GW signal. 

The IPN GRB error boxes may differ in shape and size, with areas ranging from fractions 
to hundreds or even thousands of square degrees. The short GRBs we chose for analysis have 
either small or very narrow and elongated error boxes which will make the gravitational waves 
search much easier and less computing intensive. These error boxes are tiled by a set of points 
spaced by 1.8° in each direction. An empirical upper limit of 100 deg^ square degrees was chosen 
for the maximum AA area of the error box, necessitating a maximum of about 30 independent 
sky points to search the error box. Searching over more than 30 points requires significant 
computational cost and much of the sensitivity improvement for the GRB triggered search over 
the all sky searches that have been performed would be lost. 

The short GRBs for which GW data from only the HI and H2 detectors is available are a 
special case. Since these two detectors are co-located and co-alligned, they would observe any 
gravitational wave signal at the same time, irrespective of the location of the signal. Hence, 
all the search points are degenerate, i.e. using just these two detectors would not allow us any 
spatial sensitivity since there is no time delay between these and triangulation is impossible. A 
limited size error box is not a requirement for these bursts any more and any short burst, no 
matter how extended the error box, as long as it has available data from HI and H2, will be 
analysed. Furthermore, although all-time all-sky searches for GW during S5/VSR1 have been 
published \33\ I34| . these searches did not make use of the H1-H2 data in the way we will do. 
This was because, as the detectors are co-located, they share many common sources of noise. 
Consequently, the usual method of estimating background by introducing an artificial time shift 
between the detector data is not applicable. This renders an all time search difficult as we 
have no accurate way of measuring the noise background. For the previous searches only the 
few loudest H1-H2 coincident events were considered based on no background estimations and 
solely on the coincidence test. For a GRB search, we can make use of data away from the time 
of the GRB (without time shifting) to estimate the background, therefore providing us with a 
significant increase in sensitivity. 

Depending on detector data availability and error box size we divide the GRB sample to be 
analysed into two groups: 14 short bursts with error boxes smaller than 100 deg^ and available 
data from at least two sites and 6 short bursts with an arbitrarily sized error box area that have 
available data from H1H2 only. Six bursts have large error boxes and will be considered only for 
an archival look-up in the S5/VSR1 all-sky all-time data and three other bursts have already 
been analysed and published previously. This data is summarised in Table [T] 



1. Short IPN GRBs with A A < 100 deg^ and data from two or more GW detector sites that 

will be analysed 
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2. Short IPN GRBs with data from HlH2-only that will be analysed 
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3. Short IPN GRBs with AA > 100 deg^ and data from two or more GW detector sites for 

archival data look-up only 
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4. Short IPN GRBs that have already been analysed and published 
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Table 1. The short S5/VSR1 IPN GRB sample - 14 with data from multiple non-HlH2-only 
GW detectors and well localised bursts (error box area AA <100 deg2); 6 HlH2-only poorly 
localised bursts; 6 multiple GW detectors for poorly localised bursts; 3 bursts previously analysed 
and published. At represents the time of arrival error. The IPN satellites that observed the 
bursts: S - Suzaku, Sw - Swift, I - INTEGRAL, M - MESSENGER, MO - Mars Odyssey, K - 
Konus-WIND, H - HESSI (RHESSI). 



4. The proposed search 

For the IPN short GRBs, the data streams from the operational detectors wiU be combined 
coherently and searched using the methods described in fC^. The search for compact binary 
coalescing signals is done using match filtering |36J by correlating the detector data against 
theoretical waveforms that replicate the signal for a broad interval of binary parameters. In 
order to cope with the affects of non-stationary, transient noise "glitches" in the GW detectors' 
data, the pipeline uses a number of signal consistency tests, including the null stream, amplitude 
consistency and several tests [371 [121 [3H] • The pipeline can do the search on a single sky point 
as well as on multiple sky points. Simulated GW signals injected in the data streams are used 
to assess the sensitivity of the search. The two most important changes specific to the analysis 
of the IPN GRBs are the way we will perform the search over multiple sky points and the 
generation of simulations over the GRB error boxes. 

We will generate a grid of search points to pave each GRB's error box in order to increase the 
chances of finding a signal over the entirety of the box. A fixed grid spacing (distance between 
adjacent grid points) of 1.8 degrees will be used for the search, motivated by the GW detectors' 
power of resolving the sky location, described above. Each error box is a 3-0" region but we will 
search and assign equal detection probability to each search point. 

Aside from generating a search grid, we will also generate random positions in right ascension 
and declination within an IPN error box for simulations. To generate these positions, the error 
box is first paved with grid points with a denser 0.2 degree grid spacing. This is a finer spacing 
than the one which is used for the search to ensure efficient coverage of the error box by the 
simulated positions. It will also provide a verification that the 1.8 degree spacing of search points 
is adequate. Random positions are generated within a small square bin centred at each of these 
grid points and whose sides have lengths equal to the grid spacing. The relative number (or 
density) of positions generated for each grid point is weighted according to the estimated source 
position probability distribution. The probability distribution that will be used in the case of a 
single 3-a IPN annulus is a one-dimensional Gaussian distribution centred at the central radius 
of the IPN annulus, and which has a sigma of 1/3 the width of the given annulus half- width. 
For error boxes which are formed by the intersection of two 3-a IPN annuli, the probability 
distribution will be a two-dimensional Gaussian, with distances measured from the two central 
radii of the two annuli. This assures that there are proportionally more simulations for those 
positions with larger probabilities of having a signal. Then we will draw random locations for 
injections from this list of simulation points. 

5. Discussion 

We have presented the methodology for a search for gravitational waves around the times of 
short GRBs, detected by IPN during S5 and VSRl. This search has all the needed tools and 
will commence in a short time. The work in this paper was presented in the form of a poster 
at the 10th Amaldi Conference for Gravitational Waves organised in Cardiff, UK in July 2011. 
These types of searches are very promising for the future detection of gravitational waves and 
combined with the prospect of detecting other electromagnetic counterparts from GRBs, e.g. 
radio or optical pre- or afterglows (summarised in |39| HO] ) , may open the doors for true multi- 
messenger astronomy with gravitational waves. 
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